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The theory developed by Gruner and Welsch [Phys. Rev. 
A 54, 1661 (1996)] for calculating the ID input-output rela- 
tions of the quantized electromagnetic field at dispersing and 
absorbing dielectric (multilayer) plates is generalized to the 
three-dimensional case. First a general recipe for the deriva- 
tion of the reflection and transmission coefflcients at an ar- 
bitrary body that separates two half spaces from each other 
is presented. The general theory is then applied to the case 
of planar multilayer structures, for which the Green tensor is 
well-known. 



I. INTRODUCTION 

Dielectric bodies such as beam splitters are vital ingre- 
dients in quantum-optical experiments. Their properties 
are of utmost importance for experiments with nonclas- 
sical light. For example, losses due to absorption are 
known to degrade entanglement [Q. The action on light 
of macroscopic bodies is commonly described in terms 
of input-output relations. Recently, the ID input-output 
relations of the quantized electromagnetic field at a mul- 
tilayer dielectric plate |^ have been used for constructing 
CP maps for light interacting with a dispersing and ab- 
sorbing beam splitter . 

In the present article we generalize the formalism given 
in Q and derive the 3D input-output relations of the 
quantized electromagnetic field at a dispersing and ab- 
sorbing multilayer dielectric plate. The theory is es- 
sentially based on the quantization of the macroscopic 
Maxwell field in dispersing and absorbing dielectrics (for 
a review on this topic, see Ref. [Q) where absorption is 
taken into account by introducing (for each frequency uj) 
a bosonic operator noise current density which in turn 
is related to some fundamental (collective) bosonic ex- 
citations of the electromagnetic field and the dielectric 
matter as 



^e"(r,c.)f(r,c.), 

TT 



(1) 



with e"{r,Lu) being the imaginary part of the complex 
permittivity function e(r,<j;)]. Solving the corresponding 
inhomogeneous Helmholtz equation for the electric-field 
operator in terms of the classical Green tensor G(r, s, cu), 
one obtains 



as the representation of the (Schrodinger-) operator of 
the electric-field strength in the frequency domain. 

Let us assume that the (classical) Green tensor of the 
electromagnetic field in the presence of the dielectric 
body under study is known. The input-output relations 
can then be derived by identifying the contributions to 
the electromagnetic field in the two half spaces with the 
input and output fields and the (noise) field produced by 
the absorbing body. In Section |l| we give the general 
recipe, which we specify in Section |lll| for a multilayer 
dielectric plate. 



II. IDENTIFICATION OF FIELD PARTS 

Let us assume that the space can be subdivided into 
three separated regions I, II, and III, where the region II 
represents a dielectric body. Although we shall restrict 
our attention to an infinitely extended plate of thickness 
L, the general results of this section remain also valid 
for other shapes of the body. Note that the assumption 
of infinite extension of the plate can be justified by the 
observation that, typically, the cross-sectional area of an 
impinging light beam does not exceed the surface area 
of the plate. For the sake of definiteness, let us denote 
the regions outside the body with I {z<~L/2) and III 
iz>L/2) (see Fig. |l|). 
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Ei{r, uj) = iufiQ / d'^s G(r, s, w) j(s, w) 



(2) 



-L/2 L/2 z 

FIG. 1. Scheme depicting the dielectric slab located in 
region II (—L/2 < z < L/2) surrounded by vacuum in regions 
I and III. 

The electric-field operator (|^) in region I is given by 
E(-'\r,uj)=iujfio f d'sG«(r,s,c^)j(s,c^), (3) 
where the Green tensor can be split into four parts 
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G^^\r,s,uj) = G^^"\r,s, uj) (s G I) 

+ G(")(r,s,c^) (sel) 

+ G(i2)(r,s,c^) (sell) 

+ G(i3)(r,s,c^) (selll) 



(4) 



Here, G*^^"-* (r, s, is the solution of the inhomogeneous 
Helmhohz equation (for the case where region I extends 
the whole space to infinity), and the G^^*-* (r, s, lu) are so- 
lutions of the homogeneous Helmholtz equation which 
insure the correct boundary conditions at the surfaces of 
discontinuity. The electric-field operator E'^^^'(r, w) in 
region III is given accordingly. 

Let us consider the left surface. The identification of 
the contributions of input and output fields in Eq. (|^) 
[together with Eq. (^] is unique. The term determined 
by the free Green tensor G*^^"^ (r, s, cj) represents the in- 
put field. The term determined by G*^^^^ (r, s, tj) is the 
contribution of the field which is reflected at the sur- 
face, whereas the term determined by G^^'^^ (r, s, w) de- 
scribes the transmitted field through the dielectric body 
from the right (region III). The reflected and transmit- 
ted fields can then be combined to the output field. It 
is worth noting that, since the Green tensor contains 
the full information about the allowed electromagnetic- 
field structure, the output field may also contain surface- 
guided waves that propagate along the surface. In partic- 
ular, surface-guided waves can appear if the permittivity 
of the material satisfies the conditions e'(s,cj)<0 and 
|e'(s,c.)|«|e"(s,a;)|. 

From the above, we can subdivide the field at the left 
surface (notation tc\z=-l/2 = rj) into the input field and 
the output field as 



Ei(r7,Lj) 



iir^ri,.) 



E;°"\rl,.), (5) 



with 



E 



(in) 



(r_L, w) = iujfio 



,w)j(s,w). 



(6) 



E_|°''*\r^,tj) = iw^o^c;^sG("Hr^,s,tj)j(s,w) 
Ju 

+ ioj^iQ / SsQ^^^\vl,s,oj)](s,u}). (7) 
Jill 



Analogously, for the field at the right surface (notation 

V2 = rl) 



rU=L/2 = rl) we have 



E„i(r+,^)=E;:7(r+,c.)+E::7(r|,c.), (8) 



(out) 



-III 



III 



with 



t,^'^\vl,L,)=^uJ^,o I d3sG(30)(r+ s,c.)j(s,c.), 
Jiii 



E (r+ , c.) = ^0J^,, jSs G^''^ (r+ , s, c.) j(s, c.) 

+ iuj^o / c?^sG'^^^^(r+,s,w)j(s,w) 
^11 

, / d3sG(33)(r+,s,c^)j(s,c^). (10) 
Jiii 



The arrows indicate the propagation direction of the re- 
spective field parts to larger {-^) or smaller (^) values of 
z. Recall that the output fields may also contain surface- 
guided waves whose extensions in the directions indicated 
by the arrows are small. 

The first terms in Eqs. (|^) and (|l^ are the reflected 
fields in the respective regions, the last terms are the 
transmitted fields from the opposite sides of the body, 
whereas the second terms arise from sources inside the 
body. In order to derive the input-output relations, we 
have to rewrite Eqs. (Q) and ( p^ in terms of the input 
fields at both sides of the body and noise sources from 
inside the body. This is always possible in a linear the- 
ory considered here, because the superposition principle 
holds. Thus, we can (formally) write 

E|°"'\r:,c.)= f d2s:Ri(rI,sI,c.)E;'"^(s:,c^) 



d^s+Ti,iii(r^,s+ cj)Ejjj (s+,cj) 



-Gi,i(r_L,w) + Gi,2(r^,t^), 



(11) 



E 



(9) 



;r(rl,-)= J d2slR„i(ri,s+,c.)E;;;\s+,c.) 

(fisl Tiii,i(r+, s;[, t^)^|'"' (si, w) 

+ GiiLi(r+,cj) + Giii.2(r+,cc;), (12) 

with reflection coefficients R and transmission coefficients 
T (actually second-rank tensors). The integrations run 
over the respective surfaces of the body. For example, 
the field at a point is created by the input field from 
the left hitting the surface at points leading to the 
first term in Eq. (pT|), and by the input field from the 
right hitting the surface at points thus producing the 
second term in Eq. (11). As already mentioned, surface- 
guided modes are generically included in the formalism. 
The operators G(ijn).(i,2) are related to the left and right 
propagating noise excitations inside the body. 

The reflection and the transmission coefficients in gen- 
eral depend on the observation point where the field is 
computed. For example, Ri(rJ,sJ,a;) obeys the Fred- 
holm integral equation of the first kind 

/ d2slRi(rI,s;,c.)G(i°)(sI,s,c.) - G^"' (rj, s, c.), 
J (13) 

which has to be inverted to find Ri(r2,s2,w). Post- 
multiplying Eq. (|l|) by [G^^°\s , Xj_,cli)] ^, integrating 
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over s, using the relation 



and renaming x_|_ as Sj_, we obtain the reflection coeffi- 
cient as 



Ri(r^,s^,cj) 



(15) 



Similarly, the transmission coefficient reads 
Ti.in(r2,s|,tj) 

d3sG(i3)(r-,s,c^) [g(30)(s,s+ c^) 



III 



(16) 



Analogously, the reflection and transmission coefficients 
for region III are derived to be 



Riii(r+,s+,cj) 



d3sG(33)(r+,s,c^) [g(30)(s,s+,c^) 



III 



III. I 



(r+,s^,w) 



d3sG(3i)(r+,s,L^) [g(3")(s,sX,c^) 



(17) 



(18) 



Note, that the reflection and transmission coefficients are 
polarization-dependent. The remaining task is to invert 
the free Green tensors G^^°''(s, s^, w) and G^"^"^ (s, s|[|, w). 
This can be done by expanding the inverse Green tensors 
in terms of a complete set of orthogonal solutions of the 
Helmholtz equation, i.e. the TE and TM vector poten- 
tials, and using their respective orthogonality relations. 



III. MULTILAYER DIELECTRIC PLATES 

Beam splitters and related (passive) optical elements 
typically consist of layers with different dielectric proper- 
ties (for example, anti- reflection coatings). Let us apply 
the formalism developed in Section |^ to the calculation 
of the input-output relations of light at a multilayer di- 
electric plate. To do so, we have to specify the Green 
tensor. The Green tensor for planar multilayers can be 
found in and is presented in the Appendix. All 

terms contributing to the Green tensor can be given in 
terms of TE and TM vector potentials thus making the 
distinction between different polarizations very easy. 

Suppose the plate consists of N dielectric layers (hence, 
region II is subdivided into N subregions Hi, i = l . . . N). 
Then, the input fields are given by Eq s. (P ) and (^), and 
the output fields are given by (|^) and (|10[) with the iden- 
tifications 



/ d3sG(i2)(r2,s,c^) = V / d3sG(i2*)(rl,s,c^), (19) 
[ d3sG(32)(r+,s,L^) = V / d3sG(32.)(j.+ . 

^11 j^;^ Jlli 



,s,c^), (20) 



leading to the reflection and transmission coefficients ac- 
cording to Eqs. ( p^ ) - (p^, with the Green tensor being 
specified in the Appendix. Note, that for planar multi- 
layers different polarizations do not mix, as it can be seen 
from the structure of the Green tensor, which contains 
only dyadic products of vector functions of one type. This 
behaviour is typical of planar layers. 

From the Green tensor as given in the Appendix, it is 
not difficult to identify the noise terms in Eqs. ( ITTl ) and 
( p^ ) . For propagating waves it is the subdivision into left 
and right moving waves in the slab. We derive 



Gi^i{r^, uj) 



An 



N f foo oo ( 



2 - Spr. 



Al?*Me„;,(r^,/ii)M = „;,(s, 



+ A]v'*Ne„^(rl,/ii)Ne„^(s,-/i2.) j(s,c.) , (21) 



Gi.2(r^,t^) 



47r 



j: d^s dxj:\ 



2 — Son 

\h2i 



Sl,fMe„^(r^,/ii)Me„^(s,/i2.) 



B^2^N.„;,(rI,/ii)N.„,(s,/i2,) Ms,u:)\, (22) 



Gnu(rl,-) = -^E/d^s( dXY^ 



N 



^Ili 



n=0 



2 — 5on 
\h-2t 



X [c|2»M.„^(r+,-/i3)M.„,(s,-/i2.) 
+ CrN = „A(r+,-/i3)Ne„^(s,-/i2,)l j(s,c.) \, (23) 



Giii,2(r|,t^) = - 



47r 



j: d^s dxj:\ 

^=lii. n=0 
ill 



2 — Spn 

Xh2i 



D^fMe„;,(r+,-/i3)Me„^(s,/l2 



+ i?^^^N.„^(r+,-/.3)N = „;,(s,/i2,)J j(s,c.)|. (24) 

Let us briefly compare the 3D input-output relations 
derived here with ID input-output relations given in 
Because of translational invariance along the {x,y)- 
directions and the assumed horizontal incident (propa- 
gation along the z-axis), the reflection and transmission 
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coefficients in [g| do not depend on spatial coordinates 
and reduce to scalar functions (of frequency and material 
parameters). Translational invariance along the {x,y)- 
directions is also assumed in the 3D treatment of the 
problem, but here the reflection and transmission coeffi- 
cients must be regarded as space-dependent second-rank 
tensors in general. Further, the output fields do not only 
contain volume waves but also surface-guided waves. 



IV. CONCLUSIONS 

In this article we have developed a general concept for 
deriving the 3D input-output relations of optical fields 
at dispersing and absorbing bodies, and we have applied 
it to a multilayer dielectric plate. Since the theory is 
based on the quantization of the electromagnetic field in 
absorbing media which relies on a source-quantity repre- 
sentation of the field in terms of the (classical) Green ten- 
sor, the task reduces to the determination of the Green 
tensor and its respective contributions to the input and 
output fields. In this way, the output fields can be re- 
lated to the input fields via transmission and reflection 
and to some noise excitation associated with material ab- 
sorption. Both volume waves and surface-guided waves 
are included in the input-output relations. 



APPENDIX A: GREEN TENSOR FOR A 
DIELECTRIC LAYER 

We briefly repeat some basic formulas for the Green 
tensor in the spectral representation (for details, see [^). 
The cylindrical TE and TM vector wave functions for 
even and odd waves which are frequently used throughout 
are defined as 



M = „A(r,/i) 



nJn{Xr) f sin\ 
r \cos / 



dJnjXr) 
dr 



cos 
sin 



nip 



Ahz 



(Al) 



Standard representations of the free Green tensor 
Q^^^\r,s,u>) can be found, e.g., in The scatter- 

ing Green tensor is given in ||] in the compact form of 



■roc oo , 



2-S< 



On 



(l-,5/3)Me„^(r,/j/) 

(l-5,l)4|Me„;,(s,-/l,) 

+ (l-M-B^'Me„^(s,/i,) 

X (l-,5,i)Aj;^N.„^(s,-/^3) 

+ {1~Ss3)bI;-N yj^is, hs) 
+ (l-(5/i)M.„;,(r,-/i/) 

(l-<5,i)C{;M.„;,(s,-/i,) 
+ (l-(5,3)i?{;Me„^(s,/i,) 



(A4) 



where the scattering coefficients Bj^j j^, 



and D'^j are determined by the boundary conditions 
at the interfaces between the layers. The notation used 
here is such that the indices of the Jield and source points 
cover the range (1, 2i, 3), i = 1 . . . N, distinguishing be- 
tween region I (index 1), region III (index 3), and N 
layers Hi (indices 2i) that build up region II. 



Ne„^(r,/i) = - 



dJ„(Ar) /cos\ 



ihJn{Xr) ( s 



\cos 



dr \ sin 



nipe^l, + }?Jn{Xr) 



( cos 



nil) 



where h = \/k'^ — A^, — €{lo){ijj / cY , with e(u;) being the 
complex permittivity in the respective region of space. 

The Green tensor at source point s and jield point r 
can always be decomposed into a sum of an (unbounded) 
free Green tensor G'--''"' (r, s, w) and the scattering Green 
tensor G''^^-'(r, s, w) as 

G'/^) (r, s, Lo) = G(-^o) (r, s, Lo)Sf, + G^^^' (r, s, u). 
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